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components  of  the  drag  force  on  the  particles 
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P 

R 

Re 


w 

XfV  rZ 
X,Y,Z 
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mass  of  one  particle 
gas  viscosity 

unit  vectors  in  the  direction  of  x,y,z  axes 
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velocity  in  y  direction 
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velocity  in  z  direction 
axes  fixed  in  B 
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blade  angular  velocity 
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Subscripts 

1 

2 

g 

i 

n 

out 


before  collision 
after  collision 
gas 

inlet  to  stator 
normal  to  surface 
outlet  from  stator 
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Subscripts  (continued) 


p  particles 

t  tangent  to  surface 


X 
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X 

direction 

y 
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direction 

z 
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Superscripts 


B 

Bo 

E 

P 
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relative  to  B 

the  origin  of  the  axes  xry,z  fixed  in  B 
relative  to  E 
particle  position 

a  point  that  coincides  with  p  at  a  certain  instance 
of  time 


ABSTRACT 


The  trajectories  of  solid  particles  entrained  by  fluid  flow 
in  turbine  and  compressor  stators  were  calculated.  The  impact 
and  rebound  phenomena  of  the  solid  particles  was  investigated 
experimentally  and  then  considered  in  the  solution  of  the  equations 
of  motion  of  the  solid  particles.  The  effect  of  the  particle  mean 
diameter,  material  density,  and  initial  particle  and  gas  velocities 
on  the  dynamic  behavior  of  the  solid  particles  through  the  cascade 
channel  was  also  investigated.  In  addition,  this  study  yielded 
information  concerning  blade  erosion  damage. 


INTRODUCTION 


This  investigation  was  undertaken  to  determine  a  theoretical 
approach  to  calculate  the  dynamic  behavior  of  solid  particles  en¬ 
trained  by  a  gas  flow  in  a  two-dimensional  stationary  cascade  of 
a  turbine  or  a  compressor  blading.  By  dynamic  behavior,  it  is 
meant  the  x-y  location  of  the  particles  at  any  time  in  the  cascade 
channel,  that  is,  the  path  they  follow  and  their  velocity  history 
through  the  passage  taking  into  consideration  both  their  collision 
with  the  blade  surface  and  subsequent  rebound  from  it. 

The  equations  of  motion  of  solid  particles  moving  in  a  stream 
of  gas  were  formulated  assuming  that  the  only  force  exerted  is 
the  drag  force  which  causes  the  particles  acceleration.  In  order 
to  determine  the  solution  to  the  equations  of  motion,  information 
about  the  gas  properties  in  the  channel,  the  drag  coefficient  tor 
spherical  particles,  the  restitution  ratio,  and  the  rebound  to 
incidence  angle  ratio  were  required. 

The  properties  of  the  gas  were  calculated  at  all  points  of  a 
square  grid  by  solving  numerically  the  equations  of  motion  of  gas 
flow  in  a  cascade.  The  gas  properties  were  assumed  to  be  unchanged 
due  to  the  presence  of  solid  particles. 

Drag  coefficient  formulas  for  spherical  particles  were  determined 
for  different  ranges  of  Reynolds  number  based  on  particle  mean  diameter 
and  the  velocity  difference  between  the  gas  and  particle.  Extensive 
experimental  data  were  collected  to  study  the  collision  and  rebound 
phenomena  of  the  particles  with  the  blade's  wall. 

Formulas  to  determine  the  restitution  ratio,  that  is,  the  ratio 
between  particle  velocities  after  and  before  collision,  and  the 
rebound  to  collision  angle  ratio  were  required.  This  is  necessary 
because  the  points  of  collision  are  points  of  discontinuity,  and 
hence,  new  initial  conditions  are  required  to  continue  the  solution 
of  the  equations  of  motion.  These  experimentally  correlated  formulas 
are  dependent  on  both  the  blade  material  and  the  particle  material 
and  shape,  and  may  have  different  forms  foi  different  combinations 
of  these.  However,  the  order  of  magnitude  of  these  ratios  are  the 
same,  a  fact  that  ensures  the  reliability  of  the  results  of  this 
investigation.  Once  the  operating  conditions  of  an  engine  are  known, 
these  formulas  can  be  accurately  determined  from  simple  simulated 
test  conditions.  The  equations  of  motion  of  the  particles,  as 
derived,  are _ independent  of  the  particle  concentration,  however 
they  are  valid  for  small  particle  concentrations  which  is  usually 
true  in  gas  turbine  practical  applications. 

The  effect  of  flow  parameters  such  as  particle  mean  diameter, 
particle  material  density,  initial  ncndlmensionul  particle  velocity, 
and  initial  gas  velocity  on  the  particle  trajectories  and  velocities 
are  discussed. 
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The  calculation  of  the  energy  lost  due  to  impacts  and  the 
frequency  of  impacts  at  each  location  on  the  blade  surface  are 
important  factors  in  creating  an  erosion  model  of  the  blades. 

This  information  can  be  determined  using  the  results  of  this  study , 
and  correlated  with  experimental  erosion  test  data  can  be  used  as 
a  helpful  means  to  design  turbomachine  blades  for  optimum  erosion 
as  well  as  aerodynamic  characteristics. 

The  study  of  particle  trajectories  in  stationary  flow  fields 
in  pipes  and  nozzles  have  been  reported  in  many  articles,  such  as 
References  1,  2  and  3.  However,  none  of  these  investigations  take 
into  consideration  the  collision  of  particles  with  the  wall  or  the 
cascade  type  flow.  In  Reference  4  an  experimental  study  was  made 
by  the  authors  to  determine  the  dynamic  behavior  of  solid  particles 
in  stationary  cascades.  The  patterns  of  their  trajectories  and 
velocities  were  the  same  as  predicted  by  the  present  study.  The 
locations  of  the  blade  surface  subjected  to  severe  erosion  as 
determined  from  this  study,  agrees  with  the  physical  observations 
of  erosion  tests  made  on  cascades  and  real  engine  stators.  The 
study  may  have  an  application  for  any  gas  turbine  operating  with 
particulated  gas  flow  or  any  other  type  of  two-dimensional  particle 
suspension. 


Equations  of  Motion  of  Solid  Particles 
Entrained  by  the  Gas  Flow 


Drag  force  on  a  spherical  particle  moving  in  a  stream  of  gas; 

Consider  the  motion  of  a  spherical  particle  in  a  steady,  very 
slowly  moving  incompressible  gas  flow  (Figure  a) .  The  inertia  forces 
can  be  neglected  compared  to  the  viscous  forces  and,  for  no  body 
forces,  the  Navier  Stokes  equations  become, 

grad  P  --  u  V2  (V  )  (1) 

y  y 

The  continuity  equation  of  the  gas  may  be  written  as 


div  (V  )  =0  (2) 

Taking  the  divergence  of  both  sides  of  Equation  (1)  and  using 
Equation  (2)  gives 

div  grad  P  =  V  •  V (P)  =  V2  (P)  =0  (3) 
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The  boundary  conditions  for  solving  Equation  (3)  are: 


at  the  particle  surface 


0 


at  a  distance  very  far  from  the  particle 


(4) 


(5) 


Figure  a 


From  the  physical  aspects  of  the  problem,,  and  using  the  boundary 
conditions,  a  solution  of  Equation  (3)  (Reference  5)  gives  the  pressure 
on  the  surface  of  the  sphere  as : 


P 


(6) 


The  velocity  components  of  the  gas  Ug,  v„  and  wg  in  the  x,  y  and 
z  directions  can  be  determined  from  Equation  \2)  (Reference  5).  These 
velocity  components  are: 
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_  d  x 

1  _£ _ 

8  3 
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_  d  xy  d 
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y  y»  r  4r 


d  x  z  d 
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where 


r  =  y2  +  *2 


The  shear  force  can  be  determined  at  every  point  on  the  surface 
of  the  sphere  by  using  the  partial  derivatives  of  Equations  (7)  at 
the  wall.  It  may  be  shown  that  the  resultant  of  the  pressure  dif¬ 
ferences  relative  to  the  undisturbed  flow  pressure  and  the  friction 
forces  has  the  same  magnitude  at  all  the  points  and  equals 

(3  u  |V  |)  /  d 

g  1  g„  p 

Hence,  the  total  drag  force  acting  on  the  sphere  can  be  determined 
by  multiplying  the  given  resultant  force  per  unit  area  by  the  area 
of  the  surface  of  the  sphere.  This  gives 

B  -  3  ,  ,g  dp  Vg^  (9) 

If  the  particles  are  moving  with  a  velocity  Vp  in  a  stream  of 
gas  moving  with  a  velocity  Vg,  the  whole  resistance  force  on  the 
sphere  may  be  written  as,  3 

e  "  3  *  dp  (Vg  -  y  (10) 

For  higher  Reynolds  numbers.  Equation  (10)  can  not  be  considered 
valid,  a  correction  for  Equation  (10)  was  then  required  to  allow  its 
use  for  higher  Reynolds  numbers.  A  correction  factor  which  is  a 
function  of  the  Reynolds  number  was  introduced  and  has  different 
forms  depending  on  the  range  of  the  Reynolds  number.  This  correction 
factor  was  estimated  from  experimental  results  as  discussed  later. 

The  corrected  Equation  (10)  has  the  form: 


3  ■  dp  (Vg  -  V  9<Re> 


where 


g(Ke)  =  CD  H 
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D  =  F.n.  = 
1  1 


F  n,  +  F  n.  +  F  n, 
x  1  y  2  z  3 


V  =  v  n.  =  u  n.  +  v  n~  +  w  n. 

g  gj_  i  gi  g  2  g  3 

=  x  n,  +  y  n.  +  z  ii_ 
g  1  y  2  g  3 


V  =  v  n.  =  u  nn  +  v  n_  +  w  n_ 

P  Pi  i  Pi  P2  p  3 


=  x  n,  +  y  n_  +  z  n_ 
p  1  p  2  p  3 


(12) 


In  Equation  (12),  n^,  ri2  and  n3  are  unit  vectors  in  a  frame  B 
in  the  direction  of  the  mutually  perpendicular  axes  x,  y  and  z  re¬ 
spectively.  These  axes  are  fixed  in  an  arbitrary  blade  of  the 
cascade  at  a  point  B0,  as  shown  in  Figure  b.  Dr  V  and  Vp  are 
measured  in  B.  CD  is  the  drag  coefficient  on  a  spherical  particle 
at  any  Reynolds  number. 


Figure  b 


Equations  of  motion  of  solid  part icles : 

In  order  to  derive  the  equations  of  motion  of  solid  particles 
entrained  by  gas  flow  in  a  cascade,  it  is  convenient  to  define  the 
following  axes  systems. 

_B  is  a_frame  fixed  in  an  arbitrary  blade  of  the  cascade  with 
n^,  n2  and  n3  a  set  of  nonparallel,  noncoplaner,  right  handed  unit 
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vectors  in  the  direction  of  the  axes  x,  y  and  z,  respectively  (see 
Figures  b  and  c) .  The  point  BQ,  the  origin  of  B,  is  taken  to  be 
the  intersection  of  the  plane  tangent  to  the  blade  row  at  the  entrance 
and  the  blade  axial  chord  at  the  mid-plane  that  divides  the  blade 
height  into  two  equal  parts.  Further,  frame  B  moves  with  an  angular 
velocity  EwB  in  a  reference  frame  E  fixed  in  space.  In  E,  N]_,  N2  and 
^3  are  unit  vectors  in  the  direction  of  the  mutally  perpendicular  set 
of  axes  X,  Y,  and  Z  fixed  at  a  point  E0,  which  is  the  intersection  of 
the  engine  axis  and  a  plane  normal  to  it  and  tangent  to  the  blade 
row  at  the  entrance. 

Consider  the  motion  of  a  solid  particle  p  in  B.  Let  b  be  the 
position  vector  of  p  relative  to  B0,  e  be  the  position  vector  of  p 
relative  to  E0,  and  c  the  position  vector  of  BQ relative  to  E0  (see 
Figures  b  and  c) . 


Figure  c 


Definitions  and  expressions  concerning  the  motion  of  a  particle 
in  the  frame  of  axes  described  (Reference  6  and  7)  are  required  to 
derive  the  equations  of  motion  of  particles. 
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Referring  to  Figure  b  we  may  write: 

p»  rj 

1)  w  is  the  angular  velocity  of  B  in  E  defined  as; 


E-B 

(0 


=  n. 


n3  nl  + 


n. 


nl  n2  +  nl 


n2  n3 


(13) 


where 


(14) 


denotes  the  ordinary  derivative  of  nj_  with  respect  to  time  t  in 
reference  frame  E. 

2)  Any  vector  b  in  B  can  be  written  as; 


b  =  b. 
i 


n 


=  b- 


nl  +  b2 


n2  +  b2 


n. 


(15) 


3)  The  angular  acceleration  of  B  in  E  is  defined  as  the  time 
derivative  in  E  of  the  angular  velocity  of  B  in  E,  i.e.; 


E-B  _  Ed(EwB) 
a  “  dt 


(16) 


4) 

<3l»  ^2 1 


For  a  vector  b  in 
and  qn ,  we  define 


B  which  is  a  function  of  n  scaler  variables 
n  first  partial  derivatives  of  B  in  B  as; 


3 

l 

i=l 


Ob. 

9q7  nif 


r  =  1,  -  n 


(17) 


If  b  is  only  a  function  of  time,  the  ordinary  derivative  of  b 
with  respect  to  time  in  B  becomes, 


3 

l 

i=l 


db.  db. 

_ i  -  _  _ i 

dt  i  “  dt 


(18) 


which  is  true  for  any  vector  b,  a  function  of  t  in  B. 

5)  Differentiating  Equation  (15)  with  respect  to  t  in  E,  we  get; 


E,r  db.  Edb. 

db  _  _ l  -  ,  _ i 

dt  "  dt  ni  °i  dt 


(19) 
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Using  Equations  (13),  (14),  (15)  and  (18)  in  Equation  (19), 
we  get; 


P  -  fJ 

db  _  dt  ,  E-B  v 
— =  — -j—  +  (0  x 

dt  dt 


(20) 


6)  From  Figure  b 


e  =  c  +  b 


(21) 


Using  Equation  (20)  in  Equation  (21),  the  velocity  of  p  in  E 
can  be  written  as; 

E-B  _  ^de 
V  dt 

Edc  Edb 

=  St  +  “dt 

-  EvD°  +  +  ESB  x  s 

at 


or 


E-p  E  -B-.  ,  B-p  .E-B  r 

Vc  =  V°  +  \r  +  u)  xb 


(22) 


E— D 

7)  From  Figure  b,  the  bsolute  acceleration  of  p  in  E  is  a  , 
defined  as  a  function  of  velocity  of  p  in  E  as; 


EiP  =  (EVP) 


(23) 


The  relative  acceleration  of  p  in  B  in  terms  of  the  velocity 
of  p  in  B  is 


BIP  =  (BvP) 


(24) 


The  absolute  acceleration  of  BQ  in  E  is 


EHBo  =  \  (EVBo) 


(25) 
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Substituting  Equation  (22)  into  Equation  (23),  we  get, 


EaP  =  ^§t(EVB°)  +  ^t(BVP)  +  (^~(E^B)  x  b  +  ESB  x 


dt 


(26) 


Using  Equations  (16),  (20)  and  (25)  into  Equation  (26),  we  get, 


ESP  =  ^B0  +  (_d_  (B^P)  +  E-B  x  B-P] 


E-B 


E-B  ..  ,Bdb  ,  E-B 


+  a  x  b  +  a)  x  (— +  to  x  b) 


(27) 


If  B0  was  taken  to  be  a  Point  p*  coinciding  with  P,  i.e., 
b  =  0,  Equation  (27)  becomes 


E-p  _  E-p* 

d  —  a 


B-p  ,  »  E-B  Bt-,P 
a*  +  2  (o  x 


(28) 


E-p* 

where  ab  is  the  acceleration  in  E  of  a  point  p*  that  coincides  with 
p  at  the  instant  under  consideration. 

The  equations  of  motion  of  a  particle  moving  in  a  gas  stream 
were  derived  using  Newton's  second  law  of  motion,  which  states  that 
"the  absolute  acceleration  of  the  particle  multiplied  by  its  mass 
equals  the  sum  of  the  forces  acting  on  the  particle,"  and  may  be 
written  as, 


mp  Eip  -  5 


(29) 


where  mp  is  the  mass  of  one  particle  given  by: 
.  3 

ttd 

m  =  — £ — 


(30) 


find, 


Substituting  Equations  (11),  (28)  and  (30)  into  Equation  (29)  we 


ESP  =  EJP*  +  b5P  +  2 (ESB  x  BvP) 


g(Re)]  (Dv9  -  V) 

P  P 


(31) 
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Equation  (31)  describes  the  motion  of  a  solid  particle  in  a 
gas  stream,  it  can  be  further  simplified  by  evaluating  the  dif¬ 
ferent  terms  of  the  equation, 

E-p*  E-B  -  .  E-B  ,E-B 

a^  =  a  xe  +  to  x  (  to  xe) 


-nr  (-  o  n, )  x  e  +  to  n,  x  {to  n,  x  [x  n.  + 
dt  1  1  1  p  1 


+  yp  +  (zp  +  R)  R3]} 


or 


E-p*  2  -  2  ,  .  - 

a^  =  -  to  y  n»  -  to  (z  +  R)  n, 
P  ^  P  ^ 


(32) 


as  , 


Also,  a^  is  the  acceleration  of  p  in  B  and  can  be  written 


T3Mr'l  ••  mm  **  mm  **  mm 

ap  =  x„  n.,  +  y  n_  +  z  n, 

pi  P  2  p  3 


(33) 


Furthermore,  the  Coriolis  acceleration  term  can  be  reduced  to 
2(EioB  x  BVP)  =  2((-  to  x  (x  n1  +  yp  **2  +  zp  n3)  ] 


or 


2  (EtoB  x  BVP)  -  2  to  z  n~  -  2  to  y  n, 

P  Z  P  * 


(34) 


Substituting  Equations  (32),  (33)  and  (34)  into  Equation  (31), 
it  reduces  to, 


x  Hi  +  (yp  -  to2  yp  +  2  to  zp)  H2 

+  [zp  -  to2  (zp  +  R)  -  2  (0  y  ]  n3 
18  u 

=  t r-f-2  g(Re)H(ug  -  +  (vg  -  vp)n2  +  (wg  -  wp)H3] 

PP  P 

(35) 
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Rearranging  Equation  (35),  letting 


C 


18 


g  (Re ) 


(36) 


Equation  (35)  may  take  the  expanded  form. 


XP  =  C  (ug  "  V 


yp  =  C  (vg  ~  V  +  w  yP  ‘  2  u  wp 


z_  =  C  (w  -  w  )  +  w  (z  +  R)  +  2  to  v 
P  9  P  P  P 


(37) 


In  Equation  (37) ,  the  y  distance  is  measured  in  a  straight  line 
and  not  a  curve  which  limits  the  validity  of  the  equations  to  the  case 
where  the  radius  R  of  the  cascade  row  is  relatively  large.  In  this 
investigation,  concern  was  focused  on  the  solution  of  the  two-dimensional 
problem  where  the  motion  of  gas  and  particles  was  considered  to  be  in 
the  x-y  plane  of  a  stationary  cascade. 


It  is  assumed  that. 


a)  =  0 


w 


P 


w  =  0 

g 


(38) 


Substituting  these  assumptions  into  Equation  (37) ,  they 
simplify  to  the  two  equations 


=  C  ( u  —  x  ) 

(39) 

g  p 

=  C  (v  -  y  ) 

g  p 

(40) 

Equations  (39)  and  (40)  are  the  equations  of  motion  of  solid 
particles  entrained  by  the  gas  flow  in  a  two-dimensional  stationary 
cascade.  C  is  defined  by  Equation  (36)  and  is  a  measure  of  the 
particle  translational  characteristic  time. 

Equations  (39)  and  (40)  are  two  ordinary  nonlinear  differential 
equations  that  can  be  solved  numerically  for  a  given  increment  of 
time,  to  give  the  x  and  y  coordinates  of  the  particle  as  well  as  its 
velocity  components  as  it  passes  through  the  cascade  channel.  In 
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order  to  solve  these  equations  of  motion,  the  following  information 
was  required.  First,  the  Reynolds  number  dependent  function  g(Re) 
has  to  be  determined  for  all  practical  ranges  of  application  of 
particulate  flow  in  turbomachines.  This  will  give  a  more  accurate 
estimate  of  the  drag  forces  exerted  on  the  particles.  Secondly,  the 
gas  flow  properties,  mainly  the  gas  velocity  components  Ug  and  Vg 
and  the  density  Pg,  have  to  be  known  at  all  points  in  the  flow  field 
for  certain  flow  conditions.  Finally,  the  restitution  ratio  and  the 
rebound  to  incidence  angle  ratio  has  to  be  estimated  for  certain 
particle  target  material  combinations,  which  allows  the  continuation 
of  tiie  solution  of  Equations  (39)  and  (40)  after  the  particle  impact 
with  the  wall  and  its  rebound  from  it.  This  required  information  is 
discussed  later  in  the  text  and  is  used  in  the  solution  of  Equations 
(39)  and  (40)  for  the  example  cases  of  a  turbine  stator  and  compressor 
guide  vanes  which  are  given  in  this  report. 


Drag  coefficient  for  spherical  particles: 

For  a  spherical  particle  moving  in  a  stream  of  gas  the  corrected 
formula  for  the  drag  force  on  the  sphere  is  given  by  Equation  (11) , 
which  is 


D  =  3  ug  u  dp  (V  -  V  )  g (Re)  (41) 

In  Equation  (41)  the  Reynolds  number  dependent  function  g(Re) 
has  different  forms  depending  on  the  range  of  the  Reynolds  number. 

It  should  have  a  value  of  one  for  Reynolds  numbers  less  than  1.0. 

Its  value  can  be  determined  theoretically  from  the  solution  of  the 
Navier  Stokes  equations  for  Reynolds  numbers  up  to  4.  The  experimental 
data  for  drag  coefficients  on  a  sphere  has  to  be  used  to  determine 
g(Re)  for  Reynolds  numbers  greater  than  4  and  within  the  practical 
limit  of  Reynolds  number  for  particulate  gas  flow  in  turbomachines. 

From  Equation  (10)  it  can  be  shown  that  the  drag  coefficient 
of  a  spherical  particle  at  small  Reynolds  numbers  up  to  1.0  can  be 
given  by 


(42) 


where  the  Reynolds  number  is  defined  as 


Re 


a  p 
JELLS 


(43) 


The  function  g(Re)  can  then  be  given  as  a  ratio  between  the 
drag  coefficient  at  any  Reynolds  number  and  the  drag  coefficient 
at  small  Reynolds  number,  it  may  take  the  form, 
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(44) 


g(Re)  =  CD  ff=CD/CD 

The  formulas  for  the  drag  coefficient  of  spherical  particles 
as  defined  in  Equation  (44)  are  given  below  for  different  ranges 
of  the  Reynolds  number. 


CD  =  jg  0  <  Re  <  1.0  (45) 


For  the  Reynolds  number  range  1  <  Re  <  4.0,  a  solution  of  the 
Navier-Stokes  equations  of  motion  of  a  soil'd  particle  using  a  higher 
order  term  (Reference  8)  gives  the  drag  coefficient  as, 


CD  " 


24 

Re 


(1  + 


—  Re) 
16  1 


1.0  <  Re  <  4.0 


(46) 


For  Reynolds  numbers  up  to  2000,  the  experimental  data  of 
Reference  9  for  the  drag  coefficient  of  spherical  particles  were 
fitted  and  the  drag  coefficient  was  given  by 


CD  =  21.9416  Re"0,718  +  0.3240 
4.0  <  Re  <  2000 


(47) 


For  higher  Reynolds  number  cases,  which  are  less  likely  to 
occur  in  the  problem  of  gas  particle  suspension  in  turbomachines, 
we  may  write, 

CD  =  0.4  Re  >  2000  (48) 

Detailed  discussion  of  Equations  (45)  through  (48)  may  be  found 
in  Reference  10.  These  equations,  representing  the  drag  coefficient, 
were  used  to  evaluate  the  coefficient  C  in  the  equations  of  motion 
of  the  solid  particles,  Equations  (39)  and  (40).  Figure  1  gives  the 
fitted  drag  curve  as  compared  to  the  experimental  results  obtained 
from  Reference  9  and  it  shows  very  good  agreement. 


Restitution  ratio  and  rebound  to  incidence  angle  ratio  for  solid 
particles  impinging  the  blade  ; 

In  order  to  continue  the  solution  of  the  equations  of  motion 
of  the  particles  after  the  point  of  collision  with  the  blade,  which 
is  a  point  of  discontinuity,  it  was  required  to  study  the  impact  and 
rebound  phenomena  of  the  particles  and  blade  walls „  The  value  of 
the  particle  velocity  after  rebound  and  the  angle  of  rebound  can  be 
used  as  initial  conditions  in  the  equations  of  motion  after  impact. 
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The  particles  suffer  a  loss  in  momentum  due  to  impact,  a  decrease  in 
particle  velocity  and  a  change  of  direction.  The  loss  in  particle 
momentum  due  to  collision  can  be  found  from  the  loss  in  its  velocity 
which  was  termed  restitution  ratio.  Restitution  ratio  was  defined 
as  the  ratio  between  the  particle  velocity  after  and  before  collision 

<  W  • 

The  direction  of  the  particle  after  it  collides  with  the  blade 
can  be  determined  from  the  rebound  to  incidence  angle  ratio  of  the 
particle.  This  was  defined  as  the  ratio  of  the  angle  between  the 
direction  of  the  rebound  particle  velocity  and  the  tangent  to  the 
blade  surface  at  the  point  of  impact  and  the  angle  between  the 
direction  of  incoming  particle  velocity  and  the  tangent  to  the  blade 
surface  at  the  point  of  impact  (0 2/^ 1^ •  This  information  is  enough 
to  define  the  particle  condition  after  rebound.  Since  it  was  very 
difficult  to  predict  theoretical  formulas  for  these  ratios,  it  was 
required  to  determine  them  empirically  and  examine  their  behaviors 
with  the  different  flow  parameters  and  particle  and  target  materials. 
Experiments  were  performed  in  which  the  gas  particle  flow  tunnel 
described  in  Reference  11  was  used  to  inject  particles  of  different 
sizes  in  the  flow  and  photograph  their  impact  with  a  set  of  blades 
fixed  at  the  tunnel's  test  section-. 


Particle  diameters  of  300,  1000  and  2000  microns  were  used3 
(particle  material  was  corn  cups  with  a  density  pp  =  68.7  lb/ft 
and  the  blades  were  made  of  stainless  steel).  The  analysis  of  the 
high  speed  films  taken  during  the  runs  showed  primarily,  that  the 
restitution  ratio  and  rebound  to  incidence  angle  ratio  may  be 
taken  as  a  function  of  only.  Other  parameters  seem  to  have 
little  effect  on  these  ratios.  The  results  of  these  experiments 
are  shown  in  Figures  2  through  7.  Figure  2  shows  that  the  restitu¬ 
tion  ratio  is  independent  of  the  velocity  of  particle  before  impact 
Vp  ,  decrease  in  restitution  ratio  is  shown  to  be  mainly  due  to  in¬ 


crease  in  6  ]_  as  shown  from  the  three  curves  in  Figure  2.  Figure  3 
shows  that  the  restitution  ratio  is  also  independent  of  the  place 
of  collision;  the  distance  x  meas  tred  in  the  axial  direction  for 
which  the  gas  velocity  has  different  values  at  the  suction  as  well 
as  pressure  sides  of  the  blade.  This  means  that  the  restitution 
ratio  is  not  a  function  of  the  gas  velocity  at  the  point  of  collision 
However,  the  higher  the  gas  velocity,  the  higher  will  be  Vp^  and 

hence  Vp2  will  have  a  higher  value,  tut  the  ratio  Vp2/Vp-L  independ 

ent  of  location  or  gas  velocity  at  impact,  again  6^  has  the  more 
prominent  effect,  The  data  points  :Lt  the  plot  of  Figures  2  and  3 
are  for  particle  diameters  300,  1000  and  2000  microns  and  they  are 
uniformly  distributed  along  the  curves  which  means  that  the 
restitution  ratio  is  independent  of  thv  particle  diameter.  Similar 

results  concerning  the  independence  of  Vn  /Vn,  on  d  were  shown  in 

^  2  ^1  P 

Reference  12  for  particle  diameters  fro.m  13  to  300  microns. 


In  Figure  4  the  ratio  of  the  pai  tide  velocity  tangent  to  the 
surface  after  and  before  collision  (Vp^/Vp  )  were  plotted  versus 


14 


3 2  and  may  be  considered  constant  for  all  particle  diameters.  The 
dotted  line  in  Figure  4  was  used  to  represent  the  data  on  the  figure 
and  may  be  expressed  by  the  equation 

V  /V  0.95  +  0.00055  3,  (49) 

pt2  ptl  1 

The  change  in  particle  velocity  due  to  impact  is  then  mainly 
due  to  the  contribution  of  the  change  in  its  normal  velocity.  The 
change  in  the  ratio  oetween  particle  velocity  normal  to  the  blade 

surface  after  and  before  collision  (V-  /V_  )  as  a  function  of  3, 

pn2  Pnl 

was  shown  in  Figure  5  for  different  dp.  The  dotted  curve  in  Figure  5 
represents  the  data  in  the  figure  and  may  be  written  as  follows : 


=  1.0  -  0.02108  +  0.0001417  3X  (50) 

pn3 


Equations  (49)  and  (50)  are  sufficient  to  define  the  paiticle 
velocity  and  direction  after  collision  once  its  velocity  in  the  blade 
direction  and  normal  to  it  and  its  angle  of  attack  were  known  before 
the  impact.  It  may  be  convenient,  however,  to  determine  directly 
curves  that  give  the  restitution  ratio  and  the  rebound  to  incidence 
angle  ratio  as  a  function  of  the  incidence  angle  3^ . 

Figure  6  gives  as  a  function  of  3 i  and  can  be  represented 

by  the  equation 


V 


V 


.  Cot 

B1 


-1 


[( 


t2 


nl, 


V 


V 


tl 


n2 


Cot  31) 


(51) 


32/31  has  a  minimum  value  at  3i  between  55  and  60°. 

In  Figure  7  the  restitution  ratio  Vp^/Vp^,  s  ?iven  versus  3-^ 

and  it  decreases  as  3i  increases.  The  data  in  Figure  7  may  be 
expressed  as, 


VP2 

-  , 

F 

+  Cot2 

B2 

Vpl 

\x  1 

i 

+  Cot2 

B1 

(52) 


The  general  trend  and  order  of  magnitude  of  these  results  may 
not  change  much  for  different  target  and  particle  material  (results 
for  popy  seed  particle  and  plastic  blades  had  practically  the  same 
values).  Once  the  particles  and  blade  materials  are  known,  funda¬ 
mental  experiments  can  render  the  required  information  about  the 
impact  and  rebound  phenomena.  This  information  can  then  be  used  in 
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•few* 


the  solution  of  the  equations  of  motion  of  particles  entrained  by 
gas  flow. 
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Gas  flow  properties  in  the  cascade  channel: 

The  properties  of  the  gas  flow  in  the  cascade  is  discussed 
here  very  briefly.  Reference  13  contains  a  detailed  discussion. 

For  nonviscous,  irrotational,  compressible  steady  flow  in  a 
stationary  cascade  of  a  turbomachine  the  equations  of  motion  are: 

1)  Combined  momentum  and  continuity  equations . 

a2</>  .3  (Ln  R)  3  (Ln[hP<J])  ^  3(Ln[pg])  9lJ)  _ 

_  +  - 3x .  . J  37  3y  37  "  ° 

(53) 

where 


u 


9 


1  H 

hpg  3y 


(54) 


1 

hp  3x 
y 


(55) 


2)  Combined  state  and  energy  equation. 


Boundary  conditions  for  solving  Equation  (52)  are: 


a)  =  0  along  1-4 

i|»  =  1  along  1  *  —  4 ' 

b)  Flow  along  1-2  and  3-4  are  the  same  as  the  flow 
along  1  *  - 2 '  and  3' -4'  respectively  and  the  value 
of  ij>  along  1 '  —  2  *  and  3' -4'  is  exactly  one  greater 
than  it  is  along  1-2  and  3-4. 


c) 


tan  (ei) 
(s  r) 7 


along  1-1* 


(£&) 

3*  out 


tan  (3 
(s  r) 


out 

out 


) 


along  4-4' 


(56) 


(57) 


(58) 


(59) 
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Figures  8  and  9  show  the  mesh  points  in  the  cascade  of  a  turbine 
blade  and  compressor  blades  respectively.  The  dimensions  of  the 
cascades  are  given  on  the  respective  figures.  The  computer  program 
given  in  Reference  14  was  modified  to  give  the  solution  of  Equations 
(53)  and  (56),  that  is,  the  gas  velocity  in  the  x  and  y  direction, 

Ug  and  vg,  as  well  as  the  gas  density  at  every  mesh  point  and  then 
store  the  data  on  magnetic  tapes.  Thesi  gas  properties  at  the  mesh 
point  may  be  interpolated  to  give  the  gas  condition  at  any  point  in 
the  cascade  channel.  These  values  are  needed  to  determine  the  step 
by  step  solution  of  Equations  (39)  and  (40)-.  Detailed  description 
of  the  modification  of  the  program  of  Reference  14  may  be  found  in 
Reference  15. 


Numerical  Examples 

The  dynamic  behavior  of  solid  particles  suspended  by  gas  flow 
in  the  turbine  stator  and  the  compressor  guide  vanes  shown  in  Figures 
8  and  9,  were  calculated  by  solving  the  equations  of  motion  of  the 
particles,  Equations  (38)  and  (39).  In  order  to  investigate  the 
effect  of  the  different  flow  parameters  on  the  particle  trajectories 
and  velocities,  particle  mean  diameter  was  given  the  values  1000, 

200  and  40  microns,  the  particle  material  density  was  taken  to  be 
34,  68.7  and  151  lb/ft3  and  the  initial  particle  to  gas  velocity 
varied  from  0.15  to  G.6._  Calculation  of  these  cases  was  used  to 
study  the  effect  of  dp,  p  and  Vp,/Vq.,  on  the  particle  dynamic 

behavior.  At  a  distance  far  upstream  from  the  entrance  to  the 
turbine  or  compressor  stators,  the  inlet  gas  flow  conditions  were 
gas  inlet  velocity  Vg.  =  142.64  ft/sec  in  the  axial  direction,  the 

gas  density  pa  =  0.0  76  lb/ft^  and  <-he  gas  inlet  total  cemperature 
Hi 

Ta  =  60 °F.  Gas  flow  properties,  p  ,  u„  and  v  at  every  mesh  point 

9  9  9 

were  computed  using  the  modified  program  of  Reference  14.  The  com¬ 
pressor  stator  case  requires  more  computer  time  than  the  turbine 
stator  case  in  order  for  the  solution  to  converge  for  the  same  mesh 
size. 


The  computer  program  given  by  the  authors  in  Reference  15  may 
be  used  for  the  special  case  of  two  dimensional  particle  trajectories 
to  determine  the  dynamic  behaviors  of  the  solid  particle  in  the 
stators  of  the  turbine  and  the  compressor.  The  impact  and  rebound 
characteristics  were  described  using  Equations  (51)  and’ (52).  The 
drag  coefficient  for  spherical  particles  was  represented  by  Equations 
(45)  through  (48).  The  results  of  this  investigation  are  summarized 
in  the  next  paragraph. 


Turbine  stator: 


The  turbine  blading  of  Figure  8  was  used  in  this  study.  The 
airfoil  dimensions  may  be  found  in  Reference  16.  The  particles 
were  assumed  to  enter  the  stator  uniformly  at  a  distance  one  inch 
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upstream  from  the  blade  leading  edge.  The  particle  concentration  a, 
(defined  as  the  mass  flow  rate  of  particles  divided  by  the  mass  flow 
rate  of  the  mixture) ,  specifies  the  axial  and  tangential  distance 
between  particles,  and  hence,  it  determines  the  place  and  frequency 
of  collisions,  and  consequently,  the  rate  at  which  the  blade  is  eroded. 
The  equations  of  motion  of  particles,  however,  are  independent  of  a. 

This  implies  that  a  can  be  taken  arbitrarily  to  be  any  practical  value 
without  affecting  the  reliability  of  the  conclusions  of  this  study. 

For  easier  reference,  eleven  particles  were  taken  to  enter  the 
turbine  stator  at  equal  tangential  distances  apart.  Figure  10  shows 
the  trajectories  of  the  eleven  particles  for  dp  =  1000  microns, 
while  Figure  11  shows  the  corresponding  nondimensional  velocity 
distribution  of  these  particles.  Similar  curves  for  the  trajectories 
and  nondimensional  velocities  of  eleven  particles  of  mean  diameters 
equal  200  and  40  microns  are  shown  in  Figures  12,  13,  14  and  15 
respectively.  The  first  and  the  last  particles  were  not  taken 
exactly  at  the  leading  edge  but;  slightly  shifted  in  order  to  determine 
the  direction  of  motion  of  the  particle  after  collision.  From  Figures 
10,  12  and  14,  which  show  the  particle  trajectories  for  constant  pp 
and  Vn  /V,,  ,  we  may  observe  that,  in  general,  the  particles  can  be 

*1  yi 

divided  into  three  main  categories  according  to  their  initial  place 
of  impact  (a  manner  similar  to  the  conclusions  of  the  experimental 
results  of  Reference  4. 

The  first  group  of  particles  are  the  ones  that  hit  the  upper  part 
of  the  bxade  leading  edge,  travel  laterally  in  the  nozzle,  may  or  may 
not  hit  the  opposite  blade  pressure  side,  then  move  outside  the  nozzle. 
The  s.  cona  group  consists  of  the  particles  that  hit  tne  pressure  side 
of  the  blade  then  leave  the  nozzle.  Finally,  the  last  group  of 
particles  are  the  ones  that  hit  the  lower  part  of  the  leading  edge, 
travel  laterally  and  may  or  may  not  hit  the  opposite  blade  suction 
side  or  blade  pressure  side,  then  move  away  from  the  nozzle.  It  may 
be  interesting  to  note  that,  the  fact  that  a  particle  may  or  may  not 
hit  the  opposite  blade  after  hitting  the  blade  leading  edge  and  mov¬ 
ing  laterally,  depends  on  the  pitch,  the  blade  camber  and  solidity, 
w.nich  are  geometrical  considerations.  Another  factor  is  the  mechanism 
of  the  impact  and  rebound  phenomena  for  a  certain  particle  target 
material  combination.  Factors  that  measure  the  particles  path  deviation 
from  the  gas  stream  lines  are  the  particle  diameter,  density  and  gas 
and  particle  velocities.  These  remarks  may  be  verified  by  studying 
Figures  10,  12  and  14,  where  the  throe  main  groups  of  particles  can 
be  easily  distinguished.  The  oand  width  of  the  particles  going  out 
from  the  stator  narrows  as  the  particle  diameter  decreases,  and  less 
impacts  were  observed,  particularly  with  the  blade  suction  side. 

Fror  Figures  10,  12  and  14  it  may  be  shown  that  the  areas  subjected 
to  erosion  are  the  blade  leading  edge  and  pressure  sides.  (The 
suction  side  erosi.on  may  be  due  to  particles  of  higher  diameter  re¬ 
jected  by  the  turbine  rotor  and  rebounding  back  to  hit  the  rear  part 
of  the  stator's  suction  side,  as  may  be  seen  in  Reference  17.)  The 
amount  of  erosion,  however,  is  dependent  on  the  momentum  loss  due 
to  impacts  and  the  frequency  of  impacts.  Figures  11,  13  and  15  show 
the  nondimensional  velocity  distribution  of  the  particles  given  in 
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Figures  10,  12  and  14,  respectively.  They  show  that  the  particles 
experience  loss  in  momentum  due  to  collision,  the  initial  rate  of 
increase  of  the  particle  velocity  is  higher  for  particles  of  smaller 
diameter  and  that  the  smaller  diameter  particles  attain  higher  veloci¬ 
ties.  The  particles  tend  to  reach  a  corresponding  constant  velocity 
depending  on  their  respective  relaxation  time.  The  distances  after 
which  particles  tend  to  approach  a  certain  maximum  velocity  were  of 
the  same  order  of  magnitude  (for  the  given  particles  it  is  about  2.5 
inches).  In  discussing  the  effect  of  dp,  pp  and  Vp./Vg^  on  the 

particle  dynamic  behavior  a  typical  particle  from  each  of  the  three 
groups  of  particles  were  taken  as  examples  and  each  parameter  was 
allowed  to  vary  separately  to  predict  its  influence  on  the  dynamic 
behavior  of  the  particles. 

In  order  to  show  the  effect  of  particle  mean  diameter  dp,  on 
the  dynamic  behavior,  Figures  16,  18  and  20  show  the  trajectories 
of  the_three  typical  particles  for  dp  equals  1000,  200  and  40  microns 
while  pp  and  Vp./Vg^  remain  constant,  Figures  17,  19  and  21  show  the 

nondimensional  velocity  history  of  these  particles.  Figures  16  and 
20  show  that  particles  of  smaller  diameters  experience  less  lateral 
motion  after  colliding  with  the  leading  edge  due  to  their  smaller 
inertia,  as  the  gas  stream  tends  to  force  them  to  decrease  their 
path  deviation  from  gas  stream  lines,  a  fact  that  can  also  be  observed 
from  Figure  18.  Figures  17,  19  and  21  indicate  that  particles  with 
smaller  diameters  attain  higher  velocities  in  the  same  gas  stream, 
while  starting  from  the  same  initial  velocity.  Also,  the  velocity 
drop  is  a  function  of  the  respective  angle  of  impact. 

The  effect  of  particle  material  density  pp  on  the  dynamic 

behavior  of  the  solid  particles  for  the  three  typical  trajectory 

type_was  illustrated  by  the  trajectories  of  Figures  22,  24  and  26 

for  Pp  =  34,  68.7  and  151  lb/ft3,  with  a  constant  d0  =  200  microns 

and  a  cons  tan  t  Vn./Va  =  0.3.  Figures  23,  25  and  2*7  show  the  cor- 

t*i  y 

responding  particle  nondimensional  velocity  distribution.  These 
figures  show  that  particles  with  smaller  pp  have  less  deviation 
from  the  gas  stream  than  those  with  higher  pp  and  that  they  achieve 
higher  velocities. 

The  effect  of  the  particle  initial  nondimensional  velocity  on 
the  trajectories  and  velocities  for  the  three  particle  groups  were 
shown  in  Figures  28  through  33.  Figures  28,  30  and  32  show  the 
particle  trajectories  for  vpVvg.  =  0.15,  0.3  and  0.6,  a  constant 

dp  =  200  microns  and  constant  Pp  =  68.7  lb/ft^.  Figures  29,  31  and 
33  show  the  velocities  of  these  particles  everywhere  in  the  nozzle, 
nondimensionalized  with  respect  to  the  initial  gas  velocity  of  the 
free  undisturbed  stream.  The  fig  arcs  show  that  particle  with  smaller 
initial  velocity  experience  less  deviation  from  the  gas  stream  lines 
due  to  their  smaller  inertia.  The  particles  tend  to  reach  the  same 
respective  velocity  near  the  end  of  the  nozzle  irrespective  of  their 
velocity  history.  In  the  meantime,  higher  particle  initial  velocity 
relative  to  the  gas  corresponds  to  the  case  of  smaller  gas  velocity 
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relative  to  the  particles.  Conclusions  concerning  the  effect  of 
initial  gas  velocity  would  be  that,  higher  gas  initial  velocity 
for  particles  entering  the  nozzle  with  the  same  speed,  density 
and  diameter,  will  cause  the  particles  to  follow  more  closely  the 
gas  stream  lines  and  attain  higher  velocities.  In  general,  it  may 
be  observed  that  particle  mean  diameter  has  more  influence  in  the 
particle  dynamic  behavior  than  the  particle  material  density.  On 
the  other  hand,  the  initial  particle  velocity  seems  to  have  smaller 
effect  on  the  particle  behaviors  since  they  tend  to  reach  a  constant 
velocity  irrespective  of  the  value  of  Vp./Va  • 


Compressor  guide  vanes : 


The  compressor  blading  of  Figure  9  is  given  in  detail  in 

Reference  18,  as  well  as  the  airfoil  dimensions.  This  cascade 

was  used  in  this  investigation.  Twelve  particles  were  entered  in 

the  nozzle  at  a  distance  equals  0.88  inch  upstream  from  the  blade 

leading  edge.  The  trajectories  of  these  particles  for  dp  =  1000 

microns,  pp  =  6  8.7  lb/ft^  and  Vp  /Vrf.  =  0.3  are  shown  in  Figure  34. 
r  ^i  yi 


Figures  36  and  38  show  the  trajectories  of  similar  groups  of  particles 
having  the  same  pp  and  Vp./Vg.  and  for  dp  of  200  and  40  microns, 

respectively.  Figures  35,  37  and  39  show  the  corresponding  non- 
dimensional  velocity  distribution.  The  first  and  last  particles 
in  Figures  34,  36  and  38  were  taken  slightly  above  and  below  the 
blade  leading  edge  to  uniquely  define  their  direction  after  re¬ 
bound.  Due  to  the  fact  that  compressor  blades  have  smaller  leading 
edges,  the  distance  between  the  first  two  particles  and  the  last 
two  particles  of  Figures  34,  36  and  38  were  taken  smaller  than  the 
distances  between  other  particles  in  order  to  show  the  groups  of 
particles  that  hit  the  blade  blunt  leading  edge.  Particles  are 
divided  into  four  groups,  the  group  that  initially  hit  the  upper 
part  of  the  leading  edge  and  move  laterally,  may  or  may  not  hit  the 
opposite  blade  pressure  side  then  leave  the  nozzle.  The  second 
group  are  the  particles  that  move  without  touching  the  blade  surfaces. 
Because  compressor  guide  vanes  have  smaller  blade  thickness  and 
camber  than  a  turbine  stator,  this  group  was  not  observed  in  the  case 
of  a  turbine  stator.  The  third  group  of  particles  hit  the  blade 
pressure  side  and  then  leave  the  nozzle.  The  final  group  hit  the 
lower  part  of  the  leading  edge  and  move  laterally,  may  or  may  not 
hit  the  opposit  blade  suction  side  or  blade  pressure  side,  then 
leave  the  nozzle.  It  may  be  noticed  that  similar  conclusions  to 
that  for  the  turbine  stator  case  can  be  made  for  particles  moving 
in  compressor  guide  vanes  with  one  additional  observation,  that 
the  amount  of  impact  is  smaller  for  the  case  of  compressor 
stator. 


The  effect  of  dp,  pp  and  Vp^/Vg  on  the  dynamic  behavior  of 

the  particles  moving  in  a  compressor  guide  vane  are  similar  to  the 
effects  for  the  case  of  turbine  stator.  Similar  to  the  case  of 
turbine  stator,  compressor  stator  parts  that  are  subjected  to  more 
impacts  and  hence  more  erosion,  are  the  blade  pressure  side  and 
leading  edge. 
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CONCLUSION 


For  a  given  particle  blade  material  combination,  the  restitution 
ratio  and  rebound  to  incidence  angle  ratio  may  be  determined  from 
fundamental  test  data.  The  equations  of  motion  of  particles  entrained 
by  gas  flow  in  a  turbine  or  a  compressor  stator  were  solved  using  the 
formulas  describing  the  impact  and  rebound  phenomena,  proper  drag 
coefficient  and  gas  conditions.  The  solution  of  the  particle's 
equations  of  motion  gives  the  particle  trajectories  and  velocities 
in  the  cascade  channel.  The  effect  of  particle  mean  diameter  dp, 
material  density  pp  and  initial  particle  and  gas  velocities  on 
particles  dynamic  behavior  was  discussed.  A  decrease  in  d  ,  PD  or 
Vp-/Vq.  or  increase  in  initial  gas  velocity  tends  to  stretch  the 

trajectories  in  the  axial  direction.  The  particle  mean  diameter  has 
more  effect  on  the  particles  dynamic  behavior  than  the  particle 
material  density.  The  initial  gas  and  particle  velocities  have 
smaller  influence.  The  parts  of  the  blades  subjected  to  the  maximum 
amount  of  collisions,  and  hence,  erosion  damage  are  the  blade  leading 
edge  and  rear  part  of  the  pressure  side.  For  a  turbine  stator  thr°e 
types  of  particle  trajectories  were  observed  while  for  compressor 
guide  vanes  the  same  three  groups  exist  plus  a  group  of  particles 
that  does  not  hit  the  blades  at  all.  The  results  of  this  investiga¬ 
tion,  correlated  with  experimental  erosion  tests  may  be  used  to 
design  turbomachinery  stationary  blades  for  optimum  erosion  character¬ 
istics. 
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FIG. 30  PARTICLE  TRAJECTORIES  (EFFECT  OF  Vp./Vg.) 


FIG.31  PARTICLE  NONDIME  NS  IONAL  ABSOLUTE  VELOCITY 

(  EFFECT  OF  Vp./Vg.  ) 


39 


TANGENT  I 


-15  -10  -0.5  0.0  0.5  1.0  1.5  2.0  2.5 

AXIAL  DISTANCE  -  X ,  IN. 

FIG. 34  PARTICLE  TRAJECTORIES 


COMPRESSOR  GUIDE  VANE 

dp  =  1000  microns 
pp  =  63.673  lb  /  ft3 

VP/V  03 

At  =  0.0002  seconds 


.  sr.'r:::X  *  *vfM*f***** 

<•"'* '/ rvr*-*-4  •r 


. 


]  . . 

5  rto  ^05  00  0.5  1.0  1.5  2.0  2.5 

AXIAL  DISTANCE  -X,  IN. 

35  PARTICLE  NONDIMENSIONAL  ABSOLUTE  VELOCITY 


41 


J 

< 


H  10 
Z  o 
III  | 

o 


Vp/Vg.=  0.3 

At  =  0.000133  seconds 


-1.5  -1.0  -0.5  0.0  0.5  1.0  1.5  2.0  2.5 

AXIAL  DISTANCE  -X,  IN. 

FIG.  36  PARTICLE  TRAJECTORIES 


o 

c\j 


<q 


O 


o 


COMPRESSOR  GUIDE  VANE 

dp  =  200  microns 
Pp  =  68.673  lb  /  ft3 


-1.5  -1.0  -0.5  0.0  0.5  U)  L5  ZO  2.5 

AXIAL  DISTANCE  -X,  IN. 

FIG37  PARTICLE  NONDIMENSIONAL  ABSOLUTE  VELOCITY 


42 


A/  ' 


-1.5  -1.0  -0.5  0.0  0.5  1.0  1.5  2.0  2.5 

AXIAL  DISTANCE -X,  IN. 

FIG. 38  PARTICLE  TRAJECTORIES 


FIG.39  PARTICLE  NONDIMENSIONAL  ABSOLUTE  VELOCITY 


43 


